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Abstract 

The recent WMAP data represents a milestone in cosmology and helps constrain 
cosmological parameters with unprecedented accuracy. In this work we combine the 
WMAP data with previous CMB anisotropy measurements at smaller angular scales 
to characterize the shape of the CMB anisotropy power spectrum. 

We carry out a phenomenological analysis of the data. By allowing non-physical 
shapes of the power spectrum we analyse high and low frequency experiments sep- 
arately and together. We find that WMAP dramatically constrains the power spec- 
trum up to £ ~ 700. On smaller scales, the data show discrepancies that can be 
associated with experimental systematics. If we combine all types of experiments, 
the observable features in the power spectrum are in excellent agreement with the 
WMAP cosmological parameter estimation. 

This work illustrates the advantages of a model-independent approach to under- 
standing experimental systematics that might affect CMB observations. 
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1 Introduction 



In recent years high resolution measurements of cosmic microwave background 
(CMB) anisotropies have measured the CMB power spectrum beyond the first 
acoustic peak. More recently the all-sky WMAP results represent a new step 
towards high precision cosmological parameter estimation. When combining 
these measurements with large-scale structure data and distant supernovae an 
increasingly consistent picture of the universe arises. 

In order to investigate cosmological parameters in detail the experiments 
themselves have to be thoroughly understood. CMB anisotropy experiments 
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are essentially of two kinds. High frequency experiments use bolometers and 
are usually balloon-born. The calibration, pointing and beam shape are the 
main sources of uncertainty of such instruments and tend to smear features in 
the power spectrum. Low frequency experiments use interferometers and di- 
rect^ measure the CMB in Fourier space. The calibration of such instruments 
is usually well understood and they do not suffer from beam uncertainty but 
the band-powers inferred from the observations tend to be anti-correlated. 
This can lead to increased contrast in the bands and is usually well encoded 
in the correlation matrix of the band-powers. 

This work is an extension to a previous project presented in (1). We use 
our previously obtained results as an input to a Markov Chain Monte-Carlo 
(MCMC) program to estimate the shape of the CMB anisotropy power spec- 
trum using a phenomenological curve. We use the WMAP data along with 
other CMB data separated in three subsets according to the type of experi- 
ment. 

2 Method 

In recent cosmological parameter estimations based on CMB data (see (0) 
and references therein) it has been shown that the simplest theoretical models 
based on 6 fundamental parameters are sufficient to describe the universe from 
currently available data. Common assumptions are a flat universe in which 
adiabatic initial perturbations evolve with no contribution from background 
gravity waves or massive neutrinos. 

With the recently released WMAP data a more complicated cosmological pic- 
ture can be tested (|3) including a varying spectral index, an equation of state 
of dark energy etc. It is has now become common to test a parameter space 
of 11 dimensions or more. The shape of the CMB power spectrum allowed in 
such analyses is limited by physical models. It is therefore timely to investigate 
what model-independent shape is preferred by the data in order to test the 
influence of experimental systematics which might not be fully under control. 

In this analysis, we parametrise the CMB power spectrum using a sum of five 
consecutive Gaussians 

N 

£(£ + 1)CV2tt = J2^T?exp(-(e-£ i ) 2 /2a*), N = 5. (1) 

1=1 

We then use an MCMC routine to constrain the fifteen parameters of our 
phenomenological curve; AT^, li and 0{. MCMC is an iterative process. After 
a certain burn-in time during which the Markov chains reach equilibrium, the 
parameters are sampled such that their number density is proportional to their 
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Best fit 5 gaussians - minimum chi-squared to best-fit Lambda-CDM 
spectrum from MAP anaylsis (Spergel & al. 2003) 
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Fig. 1. WMAP best-fit power spectrum fitted using a sum f(x) of five gaussians 
gi{x). The fit is discrepant at very low multipoles because gaussians cannot repro- 
duce the Sachs- Wolfe plateau and at high multipoles because there are not enough 
parameters. Between 50 < / < 1250, the fit is excellent, showing the ability of our 
phenomenological curve to fit a physically calculated power spectrum. 

posterior distribution. As explained in (4J) such a routine is sensitive to the 
prior probability assumed for the parameters. Its efficiency also depends on 
the sampling strategy. If degeneracies exist among parameters it is useful to 
sample along degeneracy directions. 

We make use of previously obtained results as an input to our MCMC rou- 
tine for both the prior and the sampling strategy. We assume our previous 
results represent a 15-dimensional Gaussian which we use a prior. For the 
sampling strategy, we diagonalise the covariance matrices obtained from our 
analysis before the WMAP data release. We then sample the parameter space 
as follows: 

P^Pf + J^CijXjPj, i, J = 1, 15, (2) 

3 

where Pi are the parameters of our Gaussians, AT^ and c^. P® are the 
previously obtained best-fit values of those parameters, are the eigenvalues 
of the covariance matrices, Cy is the matrix of eigenvectors and pj are sampled 
from a Gaussian of mean zero and unit variance. This is equivalent to following 
the axes of our Gaussian prior. We further check that the samples are within 
the previously assumed priors. 

We then run our MCMC code for three subsets of data. We first separate 
high (HF) and low frequency (LF) data as listed in table 1 and then combine 
them. We always use the WMAP data. For each sample we calculate the 
phenomenological power spectrum and return its likelihood to the MCMC 
routine. The WMAP likelihood is calculated using the method described in 
(0). For the other data the likelihood for a given phenomenological model 
is defined by -21n£ = (C P B - C B V )M BB \C V B , - C B F) where M BB , is the 
Gaussian curvature of the likelihood matrix at the peak. When available we 
use the lognormal approximation to the band-powers. We marginalize over 
the reported Gaussian distributed calibration error for each experiment and 
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Table 1 

List of experimental data used in this study. We use low and high frequency 
experiments both separately and together. 



High frequency (HF) data 


I range 


Low frequency (LF) data 


I range 


Acbar (£) 


150 - 3000 


DASI Q£) 


117 - 836 


Archeops (JJ 


15 - 350 


CBI (IV) 


400 - 1450 


Boomerang 98 (8) 


50 - 1000 


VSAE Q2) 


160 - 1400 


Maxima (9j 


73 - 1161 







we include the beam uncertainties by the analytical marginalization method 
presented in (0). 



3 Results and discussion 



We run our code until we obtain around 30, 000 samples for each run. The 
number density of each parameter is proportional to its marginalised posterior 
distribution. We summarize our resuls in table 2 along with the pre-WMAP 
results for the fifteen Gaussian parameters of our phenomenological power 
spectra. The first peak is consistenly higher when including WMAP. It is 
shifted to slightly higher multipoles. The second peak is more pronounced 
than pre-WMAP when including all data. The third peak is very high for the 
LF data. This is mainly due to the VSAE band power at £ ~ 800. At high 
multipoles, no distinct fourth or fifth peaks are detected, but the damping of 
the power is present in all cases. 

In order to get the shape of the power spectrum we divide the £ - fiK 2 space 
into pixels of size Al = 10, AT =10 and count the number of phenomeno- 
logical spectra that go through each pixel. The resulting allowed shapes are 
shown in figure 2 along with pre-WMAP results. This confirms the consis- 
tency and accuracy of the constraints obtained on the first two peaks. At 
smaller angular scales, as we had obtained before, the discrepancy remains 
between high and low frequency experiments. Even if we include correlations 
and marginalise over calibration and beam uncertainties in the data the high 
frequency data show less contrasted features in the power spectrum than low 
frequency experiments. 

In particular the third peak is higher than the second peak for low frequency 
experiments. This has important consequences in terms of the estimate of cold 
dark matter in the universe as shown in figure 3. The yellow region in figure 
3 which corresponds to the low frequency experiments intersects the orange 
regions for higher values of the CDM physical density than from the WMAP 
parameter estimation (vertical gray region) . The ratio of the first to the third 
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Table 2 

1-cr constraints on the phenomenological parameters before and after the inclusion 
of the WMAP results for three subsets of data. 



PcircirnGtGr 


All Hata 


-U WA/fAP 


n.JT del lei 


-U WMAP 

VV 1V1 Jr\ L 


1 j r Qdld 


-U WA/fAP 

VV 1V1 Jr\ L 


1X1 1 


fiQ Q+ 2 -4 


79 1 Q+0. 13 


7n 7 +4.8 
7U - 7 -3.5 


7 q Qfi+0.39 




7Q QI +0.25 


ATk 


44 - U -1.9 


/1Q 77+1-28 


4 o +3.6 


/19 7O + 2.02 


Aa o + 3.6 
4b.8_ 4 3 


,0 co+0.41 
48.5z_ Q 61 


A TU 
L\±3 


/1Q Q + 2-1 


49 7S+ - 40 


m 7 +4.3 
4U - 7 -5.0 


zii na+ - 09 
41.U9_ 63 


54 - 4 -11.2 


c/i 47+O.3O 
54.47_ 44 


A T, 




QO CQ+0.18 


ot; 7+3.1 


01 ccQ+0.00 
dl.59_ 25 


07 Q+8-7 


QO «Q + 0-H 

^•oy_ .56 


AT> r 




ib - i4 -1.04 




on 7/1 + 1-46 


on c+12.7 
^ U - b -17.3 


OQ fiR+0-36 
^• bb -1.15 


o 

Zl 


ono q+6.2 


oi A Q + 0.8 

214.3_ 7 


zU4.b_ 7 g 


01 k 0+O.2 

215.8_ 06 


one + IO.8 

206.8_ 22 


01 e t+O. 3 

216.7_ 09 


0^ 


55U_ 45 


rAO C+0.1 
04o.O_2 g 


cin^+25 
5U5_ 21 


. 19 4 +o.i 


coq+25 


coo 9+1.4 

000. 4 


h 


824+12 


831.515;° 


764+ 74 


744.5l°;2 


806±*j 


803.0lg3 


u 


H45ljg 


12191 4 


11581 242 


1241I3 


1189lg 7 


11951 2 




1474+ \f 


150ll° 


1649ljg 


17351? 


ioio_ 346 


15771 2 


cri 


93.3l 4 ; 5 2 


99.0t°; 7 


90-316.2 


99.81?; 7 


OQ r)+12.7 

88-^-12.3 


QQ q+0.5 
yy -°-0.8 


<?2 


ni.2+ 2 c 7 ; 7 


103.11°* 


7 o 9 +19.3 


61.4l?; 6 


61.9+ 3 C 6 / 7 


82.815; 2 


CT3 


po ^+20.7 


< 58.8 


136.3+ 9 C 4 ; 2 


171.4l ; 2 7 


69.8+ 12 ; 8 


< 56.8 


cr 4 


n. c. 


> 248.5 


n. c. 


248.61?- 4 


n. c. 


198.51^5 


CT5 


n. c. 


182.51?^ 


n. c. 


198.9l° 9 6 


n. c. 


6O.5IJ? 



peak from high-frequency data allows for a lower CDM density owing to the 
less marked secondary features in the power spectrum. 
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Fig. 2. Shape of the phenomenological CMB power spectrum preferred by the three 
subsets of data. The first row corresponds to pre- WMAP results which were used 
as a prior for the analysis including WMAP data which is shown on the second row. 
The first column is HF data, the second column is LF data and the third column is 
all data. The green contours are 2-a confidence levels. 

4 Conclusions 

In this work, we take advantage of the power of the WMAP data to constrain 
the shape of the CMB power spectrum. We find that when including data from 
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Fig. 3. Ratio of amplitudes of the third peak to the first and second peaks vs the 
physical density of cold dark matter. Our phenomenological results differ. However, 
when combining all data, we are consistent with the WMAP estimation (vertical 
gray region). The orange region corresponds to the degeneracy direction for physi- 
cal models. The grey horizontal region corresponds to pre- WMAP results and the 
brown region corresponds to a pre- WMAP parameter estimation. Our results for 
the HF, LF and all data are shown in the blue, yellow and green horizontal strips, 
respectively. 

both high and low frequency experiments the preferred shape of the power 
spectrum is very close to the shape of the physical 'concordance model'. Low 
frequency data show very contrasted second and third peaks. High frequency 
data prefer more shallow features. In all cases a damping of the power on small 
scales is detected. 

Such a model-independent approach sheds light on discrepancies between high 
and low frequency experiments which would be less marked in traditional 
parameter estimation techniques. Our best fit shapes are not physical but the 
main observables of the spectrum are consistent with the WMAP cosmological 
parameter estimations. 

Further investigation of the experiment-related discrepancies would be useful 
and WMAP provides a strong starting point. In combination with higher mul- 
tipole data the shape of the CMB power spectrum is now strongly constrained. 
There is no doubt that a second and third peaks have been detected. This il- 
lustrates the new era cosmology has entered: the question now is not whether 
secondary peaks are present, thereby strengthening the already compelling 
agreement with inflationary predictions but how strong those secondary fea- 
tures are. 
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